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ABSTRACT 

We have conducted a survey of intermediate Galactic latitudes using the 13-beam 
21-cm multibeam receiver of the Parkes 64-m radio telescope. The survey covered the 
region enclosed by 5° < \b\ < 15° and —100° < I < 50° with 4,702 processed pointings 
of 265 s each, for a total of 14.5 days of integration time. Thirteen 2 x 96-channel 
filterbanks provided 288 MHz of bandwidth at a centre frequency of 1374 MHz, one- 
bit sampled every 125 [is and incurring ^DM/13.4 cm™ 3 pc samples of dispersion 
smearing. The system was sensitive to slow and most millisecond pulsars in the region 
with flux densities greater than approximately 0.3-1.1 mJy Offline analysis on the 64- 
node Swinburne workstation cluster resulted in the detection of 170 pulsars of which 
69 were new discoveries. Eight of the new pulsars, by virtue of their small spin periods 
and period derivatives, may be recycled and have been reported elsewhere. The slow 
pulsars discovered are typical of those already known in the volume searched, being 
of intermediate to old age. Several pulsars experience pulse nulling and two display 
very regular drifting sub-pulses. We discuss the new discoveries and provide timing 
parameters for the 48 slow pulsars for which we have a phase-connnected solution. 
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1 INTRODUCTION 

By the late 1990s radio pulsar surveys had resulted in the 
discovery of ~700 pulsars, spawning numerous studies with 
wide ranging implications for astrophysics and physics in 
general. Despite having been first discovered over a quar- 
ter of a century earlier, pulsars with unique and interest- 
ing properties (e.g. Wolszczan & Frail 1992; Johnston et al. 
1992b; Johnston et al. 1993; Bell et al. 1995; Stappers et al. 
1996) continued to be uncovered by surveys which also 
served the purpose of providing a larger sample for statis- 
tical analyses of classes of pulsars and pulsar binaries (e.g. 
Lyne et al. 1998). 

Nearly all early surveys were conducted at low frequen- 
cies (y ~ 400 MHz) due to the steep spectrum (a ~ —1.6, 
where S a !/ a ; Lorimer et al. 1995) characteristic of mi- 
crowave radiation from pulsars and the faster sky cover- 
age afforded by the larger telescope beam at these frequen- 
cies. However, two effects that hamper the detection of cer- 
tain pulsars at low frequencies can be avoided by using a 
higher frequency. Firstly, for small Galactic latitudes the 
background of Galactic synchrotron emission comprises the 
main contribution to the system temperature at these fre- 
quencies. The spectrum of this radiation is steep (a ~ —2.6; 
Lawson et al. 1987) and at high frequencies generally repre- 
sents an insignificant contribution compared to the thermal 
receiver noise. Since they share the low Galactic z-height of 
their progenitor population, young pulsars in particular are 



selected against in low frequency surveys due to the elevated 
sky background temperature. Secondly, radiation propagat- 
ing through the interstellar medium is subject to 'scatter- 
ing' due to multi-path propagation, effectively convolving 
the light curve with an exponential of a time constant that 
scales as (Abies, Komesaroff & Hamilton 1970). Since 

the minimum detectable mean flux density in pulsar obser- 
vations is proportional to [5/(1 — S)] 1 ^ 2 where 5 is the ef- 
fective pulse duty cycle, scatter-broadening of the received 
pulses hampers the detection of pulsars at low frequencies, 
especially those with short spin period such as the interest- 
ing and important class of 'millisecond' pulsars, and (again) 
young pulsars. Moreover, by conducting a survey at high fre- 
quencies one is sensitive to pulsars with flatter spectra that 
were missed in earlier surveys. 

With the rise in availability of affordable computing 
power in the 1980s it became feasible to process surveys 
with fast sampling rates and large numbers of pointings, as 
required for large scale high frequency surveys for millisec- 
ond pulsars. Clifton et al. (1992) and Johnston et al. (1992a) 
conducted highly successful 20-cm pulsar surveys near the 
Galactic plane, discovering 86 pulsars between them, includ- 
ing a high fraction of young pulsars. However, the surveys 
did not have sufficient sensitivity at high time resolution to 
discover any millisecond pulsars. In addition, for the reasons 
mentioned above the surveys concentrated on the Galactic 
plane and hence the samples of detected pulsars were of 
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Figure 1. The multibeam tessellation unit shown with circles 
depicting the half-power points of beams. A unit is observed with 
four offset pointings, one of which is hatched in the above for 
clarity. The shape made by the 52 beams can be seamlessly self- 
tcsscllated. 



reduced value in modelling the Galactic pulsar population 
compared to larger surveys. 

In 1997 the Australia Telescope National Facility com- 
missioned a new 21-cm 13-feed multibeam receiver, primar- 
ily for HI surveys (Henning et al. 2000; Barnes et al. 2001). 
The large instantaneous sky coverage and excellent sensitiv- 
ity also makes the system a powerful pulsar survey instru- 
ment and this led to the commencement of a long-running 
deep survey of the southern Galactic plane (|6| < 5°) which 
is expected to almost double the known population (Lyne 
et al. 2000; Camilo et al. 2000). We conducted Monte Carlo 
simulations similar to those discussed by Toscano et al. 
(1998) and found that a shallower 'flanking' survey should 
discover a sizeable population of pulsars with unprecedented 
time efficiency in an area of sky not previously sampled at 
high frequencies. Based on this result we conducted such 
a survey between 1998 August and 1999 August. The sur- 
vey proved highly successful, discovering 69 pulsars includ- 
ing two pulsar binaries containing heavy CO white dwarfs, 
one of which will coalesce in less than a Hubble time with 
dramatic and unknown consequences (Edwards & Bailes 
2001b), and a further four binary and two (perhaps three) 
isolated recycled pulsars with important implications for 
theories of binary evolution (Edwards & Bailes 2001a). In 
this paper we report in detail on the observing system, anal- 
ysis procedures, sensitivity and completeness. We discuss 
detections of previously known pulsars and present the new 
sample of slow pulsars, including timing results for those 
with solutions. 



2 OBSERVATIONS AND ANALYSIS 



2.1 Hardware Configuration and Survey 
Observations 

The 64-m Parkes radio telescope was used with the 13-beam 
21-cm receiver (Staveley-Smith et al. 1996) which provides 
300 MHz of bandwidth and a system temperature of ~21 K. 
Signals from the two orthogonal polarisations of each beam 
were mixed with a local oscillator before being fed to an 
array of 26 96-channel filterbanks. Each filterbank channel 
was 3 MHz wide and the band was centred at a frequency 
of 1374 MHz. The detected signals from corresponding po- 
larisation pairs in each channel were summed and high pass 
filtered (with a time constant of ~0.9 s; Manchester et al. 
2000) before being integrated and one-bit sampled every 125 
/is. The data stream was written to magnetic tape (DLT 
7000) for offline processing, as well as being made available 
to online interference monitoring software in near-real-time 
via the computer network. With the exception of the sam- 
pling interval, the system was identical to that used for the 
Galactic plane survey (Manchester et al. 2001). 

The receiver feeds are arranged in such a way as to al- 
low coverage of the sky in a hexagonal grid, with beams 
overlapping at their approximate half-power points (7' from 
the beam centre). A group of four pointings results in the 
uniform coverage of a roughly circular shape ~ 1 degree 
in radius which in turn can be efficiently tessellated (see 
Figure [l]). The region enclosed by 5° < |6| < 15° and 
— 100° < I < 50° was covered in 4,764 265-s proposed point- 
ings, amounting to only 14.6 days of integration time. Most 
of these pointings were observed in several week-long ob- 
serving runs between August 1998 and August 1999. 



2.2 Search Analysis Procedure 

The processing of the 64-tape ~ 1.6 terabyte data set was 
performed on the Swinburne Supercluster, a network of 64 
Compaq Alpha workstations. Before searching for pulsars, 
each beam was analysed for the presence of powerful signals 
that appeared in only a few filterbank channels, a common 
type of interference signal. When such signals were present, 
samples in the culprit channels were zeroed, a process that 
does not incur too much loss of sensitivity since this varies 
as the square root of the effective bandwidth. In addition, 
broad-band periodic signals that appear in large numbers 
of beams in any given 30-minute period were detected and 
logged to a file for later reference. 

To correct for the effects of non-linearity in the dis- 
persion relation, data from the 96 filterbank channels were 
padded with 32 dummy channels in such a way as to al- 
low linear de-dispersion of the resulting 128 channels, as 
used by the Galactic plane survey collaboration. This en- 
abled the use of the fast 'tree' algorithm of Taylor (1974) 
to partially de-disperse the data into eight trial dispersion 
measures in each of sixteen sub-bands. Whilst the linearity 
of dispersion with respect to channel number would enable 
full de-dispersion (that is, 128 trial DMs with no sub-band 
divisions), in order to limit storage requirements and to al- 
low the recording of frequency-resolved pulse profiles to aid 
in suspect scrutiny, application of the algorithm was stopped 
after the production of eight DMs. 

The tree algorithm produces trial DMs up to the 'diag- 
onal' DM of 17.0 cm -3 pc, where the dispersion delay across 
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Figure 2. Sky coverage for the survey. Ellipses represent groups of four inter-meshed pointings. Shades of grey represent the density of 
coverage, from unobserved (white) to observed twice (black). 

one sub-band in units of samples is equal to the number of 
channels used to form it. It should be noted that in previ- 
ous surveys where the linear dispersion approximation was 
acceptable for the tree stage, this parameter was approxi- 
mately equal to the DM at which the smearing induced in 
each channel was one sample interval. The latter parame- 
ter is commonly quoted in conjunction with the sampling 
interval to give an indication of the time resolution avail- 
able to pulsars of various DMs in a pulsar survey. For the 
present survey this value varies from 9.4 to 17.5 cm~ 3 pc 
depending on the centre frequency of the channel, and for 
evaluation purposes one should use the geometric mean of 
13.4 cm~ 3 pc. The tree algorithm was extended to also pro- 
duce time series for 1-2 times the diagonal DM, and beyond 
this value the sample interval was doubled by summing of 
samples before re-application of the algorithm, and the pro- 
cess repeated to produce time series with 2-4, 4-8, 8-16 and 
16-32 times the diagonal DM. 

The periodicity search itself was based on that of the 
Parkes Southern Pulsar Survey (Manchester et al. 1996), 
generalized and modified to handle the large number of spu- 
rious interference signals present in the multibeam data. 
Time series were constructed at 375 trial values of disper- 
sion measure from to 562.5 cm -3 pc by summing partially 
de-dispersed sub-bands in the nearest DM with the appro- 
priate time offsets. The trial DMs were spaced in such a way 
that the effective smearing induced due to the difference be- 
tween the DM of a pulsar and the nearest trial DM was no 
more than twice that induced by the finite width of individ- 



ual filterbank channels. The time series were filtered with a 
boxcar of width 2048 ms to remove the effects of receiver 
noise and gain variations during the course of the observa- 
tion, before being Fourier transformed and detected to form 
the fluctuation power spectrum. 

For signals with frequencies lying on the boundary be- 
tween two spectral bins the result is two components of equal 
magnitude and opposite sign in the adjacent bins. To main- 
tain sensitivity to such signals we also computed the differ- 
ence of each bin and its neighbours and used half the squared 
magnitude of the results as alternative estimates of spectral 
power. For each bin the highest of the three power values 
computed was chosen for use in the final power spectrum. 
In the case of the zero-DM time series, this spectrum was 
checked for the occurrence of any signal with a frequency 
close to one earlier logged as a broad band interference sig- 
nal contemporaneous with this observation. Should such a 
signal be present, its exact extent in the spectrum was as- 
sessed and the corresponding bins zeroed in this and all 
other power spectra searched in this beam. The spectrum 
above a frequency of 1/12 Hz was then searched for signifi- 
cant spikes compared to a local mean (to compensate for the 
overall redness of the spectrum). Harmonics were summed 
and the process repeated for up to 16 harmonics to main- 
tain sensitivity to signals with short duty cycles. Significant 
signals at any level of harmonic summing were recorded and 
after all trial dispersion measures had been searched the set 
of signals was correlated into a number of candidates, each 
covering signals of similar pulse period occurring at multiple 
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Figure 3. Estimated minimum detectable mean flux density 
(Smin) as a function of pulse period for intrinsic pulse widths 
of 10° (solid lines) and 90° (dashed lines) at dispersion measures 
of 0, 10, 30, 100 and 300 cm -3 pc (in order of increasing 5 m i n for a 
given pulse period). Points represent undetected pulsars which lie 
within 10' of an observed beam, where flux density measurements 
have been published. Flux densities published without uncertain- 
ties are plotted without error bars and in such cases the relative 
uncertainty is probably around 50 per cent. 



trial DMs. The top 99 candidates in each beam were sub- 
ject to a fine search (by means of maximisation of signal to 
noise ratio, S/N) in period and dispersion measure around 
the best values found in the spectral search. Pertinent infor- 
mation including the resulting best profile, grey scale maps 
of pulse profiles as a function of time and radio frequency 
and of signal to noise as a function of period and dispersion 
measure were saved to disk. 



2.3 Suspect Scrutiny, Confirmation and Timing 
Observations 

The final stage of analysis was human viewing. The large 
number of beams and the prevalence of interference signals 
presented considerable complications to the viewing process 
due to the volume of candidates produced. In previous sur- 
veys (e.g. Manchester et al. 1996) candidates of similar pe- 
riod occurring in multiple beams contemporaneously were 
generally taken as interference signals and ignored. In the 
case of results from this survey, the plethora of interference 
signals across the spectrum resulted in the misinterpreta- 
tion of many pulsars as interference. It was found that this 
limited the applicability of this approach to the handful of 
periods that appeared more than ~ 250 times on any tape. 
All remaining candidates with signal to noise ratios greater 
than eight (of which there were several hundred thousand) 
were then scrutinized by a human viewer and promising sig- 
nals scheduled for confirmation by re-observation. 

Human viewing of all suspects with S/N > 8 was ex- 




0.1 1 10 

Flux Density at 1400 MHz 

Figure 4. Histograms depicting distribution in flux density of all 
previously known pulsars with published flux density with pro- 
cessed beams centred less than 10' away. The subset of such pul- 
sars that were undetected in the survey is represented by the 
hatched regions. 

pected to be incomplete in its selection due to the viewing 
speed necessitated by the large number of candidates to be 
assessed. This method was used as a first pass over the data, 
however after all data were processed and the 30 GB result 
set assembled, a more complete candidate analysis scheme 
was employed. A custom-written graphical software package 
allowed for visual (and numerical) identification of the dis- 
tribution of candidates in a variety of parameters. Sets of 
candidates could be trimmed by the graphical or command- 
driven selection and deletion of interference signals, and the 
remaining candidate list subjected to human scrutiny. A 
set of 'macros' were developed for the deletion of dozens 
of commonly appearing highly coherent interference periods 
as well as all signals with large relative errors in disper- 
sion measure (a characteristic of terrestrial, non-dispersed 
interference) and any candidate with a S/N less than 9 (or 
9.5 for P < 20 ms). These produced an order of magnitude 
reduction in the number of suspects to be viewed, and a cor- 
responding improvement in the accuracy and completeness 
of scrutiny. Interference mitigation procedures employed for 
each tape were also recorded as macros to allow for repeata- 
bility and quantification of any selection effects imposed. 

Those candidates confirmed by detection in a re- 
observation were added to an ongoing program of pulsar tim- 
ing of new discoveries. Observations of typically 250 s were 
made with the centre beam of the system described above, 
or more recently the 2 x 512 x 0.5-MHz filterbank to provide 
improved time resolution for short period pulsars. For most 
pulsars at least one timing observation was obtained at a fre- 
quency of 660 MHz to enable accurate measurement of the 
dispersion measure. In offline processing samples were de- 
dispersed and folded at the predicted topocentric pulse pe- 
riod. The resulting pulse profiles were fitted to a 'standard' 
profile usually produced by adding several prior observations 
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Table 1. Detected previously known pulsars 
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t — Pulsar lies outside nominal survey region 

t — Pulsar originally published with incorrect period; corrected period listed 



© 0000 RAS, MNRAS 000, 000-000 



6 R.T. Edwards, M. Bailes, W. van Straten and M.C. Britton 



Table 1 — continued 
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6.0 




48.7 


B1732-07 


0.419 


73.5 


1 7(2) 

_L. 1 


17.3 


13.3 


3.5 


78.4 


B1738-08 


2.043 


74.9 


1.4(4) 


17.0 


11.3 


4.9 


91.3 


B1740-13 


0.405 


115 


0.50(10) 


12.7 


8.2 


3.9 


24.0 


J1744-1134 


0.004 


3.1 


2.0(2) 


14.8 


9.2 


3.2 


24.8 


B1745-12 


0.394 


100.0 


2.0(3) 


14.0 


7.7 


3.8 


106.1 


B1747-46 


0.742 


20.3 


10.00 


-15.0 


-10.2 


6.5 


159.0 


B1758-03 


0.921 


117.6 


0.70(10) 


23.6 


9.3 


4.1 


30.3 


B1802+03 


0.219 


79.4 




30.4 


11.7 


6.4 


27.7 


B1804-08 


0.164 


112.8 


16(4) 


20.1 


5.6 


5.2 


260.9 


J1808+00 


0.425 


141 




28.5 


9.8 


8.1 


21.4 


J1808-0813 


0.876 


151.3 


2.00 


20.6 


5.8 


7.1 


34.9 


B1810+02 


0.794 


101.6 


0.30 


30.7 


9.7 


4.7 


14.6 


J1811+0702 


0.462 


54 




34.7 


12.1 


2.2 


21.0 


B1813-26t 


0.593 


128.1 




5.2 


-4.9 


7.8 


17.9 


B1813-36 


0.387 


94.4 


2.00 


-3.2 


-9.4 


5.4 


52.4 


J1817-3837* 


0.384 


102.8 




-5.3 


-10.4 


7.4 


57.6 


B1818-04t 


0.598 


84.4 


8.0(6) 


25.5 


4.7 


12.0 


104.6 


B1820-31 


0.284 


50.3 


2.5(6) 


2.1 


-8.3 


8.3 


72.1 


B1821+05 


0.753 


67.2 


1.7(4) 


35.0 


8.9 


3.3 


111.9 


B1822+00 


0.779 


54.4 


0.40(10) 


30.0 


5.9 


1.4 


26.6 


J1822-4209 


0.457 


72.5 


1.50 


-8.1 


-12.8 


6.3 


16.2 


B1839+09 


0.381 


49.1 


1.70(10) 


40.1 


6.3 


3.5 


99.9 


B1842+14 


0.375 


41.2 


1.5(3) 


45.6 


8.1 


5.2 


30.5 


B1845-19 


4.308 


18.3 




14.8 


-8.3 


9.3 


84.6 


B1848+12 


1.205 


71 


0.50(10) 


44.5 


5.9 


5.5 


21.3 


B1848+13 


0.346 


59.0 


1.4(3) 


45.0 


6.3 


8.1 


32.1 


J1848-1414 


0.298 


134.4 




19.9 


-5.8 


4.3 


11.2 


B1851-14 


1.147 


130.1 




20.5 


-7.2 


9.0 


24.0 


J1852-2610 


0.336 


56.8 


1.40 


9.5 


-11.9 


0.8 


36.2 


B1857-26 


0.612 


38.1 


13.0(10) 


10.3 


-13.5 


6.8 


201.6 


B1900-06 


0.432 


195.7 




28.5 


-5.7 


10.0 


13.7 


J1901-0906* 


1.782 


72.7 


3.10 


26.0 


-6.4 


4.7 


98.6 


J1904-1224 


0.751 


118.2 




23.3 


-8.5 


2.5 


14.9 


B1907-03 


0.505 


205.7 


0.80 


32.3 


-5.7 


3.6 


22.1 


B1911-04 


0.826 


89.4 


4.4(5) 


31.3 


-7.1 


4.8 


233.4 


B1917+00 


1.272 


90.7 


0.8(2) 


36.5 


-6.2 


4.6 


36.7 


B1923+04 


1.074 


101.8 




41.0 


-5.7 


9.0 


15.5 


J1929+00 


1.167 


33 




37.7 


-8.3 


2.3 


14.0 


J1938+0652 


1.122 


70 




44.4 


-7.1 


3.8 


29.2 


B1942-00 


1.046 


58.1 


0.80(10) 


38.6 


-12.3 


6.2 


21.5 


t — Pulsar lies outside nominal survey region 











t — Pulsar originally published with incorrect period; corrected period listed 
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Table 2. Undetected previously known pulsars 

Name P DM S1400 I b A pos 

(s) (cm-3pc) (mjy) (°) (°) (') 



B0923-58 


0.740 


57.7 




-81.6 


-5.6 


15.7 


J1006-6311 


0.836 


196.0 




-74.4 


-6.0 


21.6 


J1123-6651 


0.233 


111.2 


0.50 


-65.5 


-5.4 


5.4 


J1130-6807 


0.256 


148.7 




-64.5 


-6.4 


5.8 


J1137-6700 


0.556 


228.0 


1.10 


-64.2 


-5.2 


17.7 


J1143-5158 


0.676 


159.0 


0.30 


-67.6 


9.5 


3.4 


J1225-5556 


1.018 


125.8 


0.40 


-60.7 


6.7 


6.6 


J1356-5521 


0.507 


174.2 


1.50 


-47.8 


6.3 


11.8 


B1503-51 


0.841 


61.0 




-36.9 


5.5 


8.4 


J1604-7203 


0.341 


54.4 




-43.3 


-14.6 


i.l 


J1654-2713 


0.792 


92.4 


0.10 


-5.0 


10.3 


3.1 


B1659-60 


0.306 


54 




-30.2 


-11.4 


6.3 


B1700-32 


1.212 


109.6 


6.00 


-8.2 


5.4 


21.0 


J1701-3006 


0.005 


115.6 




-6.4 


7.3 


6.6 


J1732-1930 


0.484 


73.0 


0.30 


6.4 


7.6 


5.3 


B1740-03 


0.445 


30.2 


0.50 


21.6 


13.4 


5.0 


J1740-5340 


0.004 


71.9 




-21.8 


-12.0 


3.7 


B1745-56 


1.332 


58 




-23.4 


-14.3 


50.1 


B1802-07 


0.023 


186.4 


0.60 


20.8 


6.8 


7.6 


J1807+07 


0.464 


89 




35.1 


13.3 


i.l 


J1809-3547 


0.860 


193.8 




-3.5 


-7.8 


7.6 


B1820-30A 


0.005 


86.8 


0.72(2) 


2.8 


-7.9 


6.6 


B1820-30B 


0.379 


87.0 


0.07(3) 


2.8 


-7.9 


6.6 


J1821+17 


1.366 


79 




45.3 


14.2 


8.7 


B1821-24 


0.003 


119.8 


1.80 


7.8 


-5.6 


6.4 


J1822+0705 


1.363 


50 




36.0 


9.7 


4.2 


J1822+11 


1.787 


112 




39.9 


11.6 


7.4 


J1823-0154 


0.760 


135.9 


0.80 


28.1 


5.3 


28.7 


J1834+10 


1.173 


62 




40.6 


8.6 


8.5 


J1838+06 


1.122 


70 




37.5 


6.1 


1.0 


J1838+16 


1.902 


36 




46.7 


10.3 


10.2 


J1859+1526 


0.934 


97.4 




47.6 


5.2 


76.0 


J1911-1114 


0.004 


31.0 


0.70(10) 


25.1 


-9.6 


6.2 


J1933+07 


0.437 


170 




44.8 


-5.6 


10.2 


J1941+1026 


0.905 


138.9 




48.0 


-6.2 


15.9 


J1947+10 


1.111 


149 




49.0 


-7.3 


56.5 


J1950+05 


0.456 


71 




44.9 


-10.6 


5.5 



4,702 pointings. In the first observing run an error in the 
telescope control system resulted in spurious rotations of the 
receiver feed, making the sky position corresponding to each 
recorded beam (except the centre beam) indeterminate and 
reducing sensitivity by moving off (or on!) source midway 
through observations. For this reason most of the region of 
sky surveyed under these conditions were re-observed. 

Figure || shows the sky coverage achieved by the survey. 
One ellipse is plotted for each group of four inter-meshing 
pointings. Ellipses are shaded according to the number of 
observed and processed pointings in the group. Most groups 
are either medium grey (for standard once-only coverage), 
black (for those areas observed twice due to position uncer- 
tainty as discussed above) or white (for unobserved groups). 
Other shades reflect varying numbers of observed and pro- 



of high signal to noise ratio. The resulting phase offsets were 
used to produce barycentric times-of-arrival (TOAs) which 
were then used in conjunction with the TEMPO software 
package^ to fit for the relevant spin, astrometric and binary 
parameters (see e.g. Taylor & Weisberg 1989). 



3 RESULTS AND DISCUSSION 

3.1 Sky Coverage, Sensitivity and Re-Detections 

The survey was deemed complete in August 1999 after the 
observation and successful processing of 60,852 beams in 

* http:/ /pulsar. princeton.edu/tempo 
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Figure 5. Pulse profiles for new slow pulsars 
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Figure 5 - continued 



cessed beams in the group. A total of 4,465 of the proposed 
pointings were observed at least once, yielding a metric of 
completeness of 4,465/4,764 = 94 per cent. 

The sensitivity of a pulsar survey is derived from the 
radiometer equation, altered to take into account the detec- 
tion of periodic signals. It can be expressed as follows, 



ap(T Tec +T Bky ) ( 6 \ 1/2 
G(N p Avt ohs ) 



si ( 6 V 

V2 \ l-8j 



(1) 



where Smin is the minimum detectable mean flux density, a 
is the threshold S/N, /3 is a dimensionless correction factor 
for system losses, T rcc and T sky are the receiver and sky noise 
temperatures, G is the telescope gain, N p is the number of 
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Figure 6. Distribution in pulse period and period derivative of new pulsars (stars), previously known pulsars in the survey region (large 
dots) and all other known pulsars (small dots), where such values have been measured. Also plotted are solid and dashed lines of constant 

characteristic age (t c = P/2P) and inferred surface magnetic field strength (B = 3.2 X 10 19 G s -1 / 2 \J PP) respectively. 



polarisations, Av is the observing bandwidth, t b s is the in- 
tegration time, W is the effective pulse width in time units 
and P is the pulse period (Dewey et al. 1985). The effective 
pulse width (W = 5P where 5 is the observed duty cycle) is 
computed as the quadrature sum of the intrinsic pulse width 
and pulse broadening terms due to such effects as dispersion 
smearing, scatter-broadening and the finite sampling inter- 
val. 

The system characteristics of the multibeam receiver 
vary from beam to beam and we use here averages of the 
values presented at the instrument web page^], yielding a re- 
ceiver temperature of 21 K and a gain of 0.64 K Jy _1 . The 
dimensionless parameter f3 embodies the loss due to one-bit 
digitisation (\/V/2 ~ 1.25) and other system losses, which 
we treat collectively with f3 — 1.5. Assuming a typical sky 
temperature of 2 K, the calculated sensitivity as a function of 
pulse period is plotted in Figure ^ for several dispersion mea- 
sures and pulse widths. Included in the effective width calcu- 
lation are the dispersion smearing in filterbank channels and 
the sampling interval. It should be noted that the sensitivity 
derived is that available at the centre of the beam. The aver- 



t 

http: / / www.atnf.csiro.au / research / multibeam /lstavele / description 



age beam response is approximately Gaussian (Manchester 
et al. 2001) with a half-power width of 14.3', meaning that 
pulsars lying near the beam overlap points will be detected 
with half the sensitivity calculated. We calculate that the 
mean sensitivity within the half-power radius is 73 per cent 
of that at the centre of the beam, meaning that for statis- 
tical sensitivity estimates, the curves in Figure ^ should be 
raised about a tenth of a decade. The resulting basic mean 
sensitivity to slow pulsars is 0.3-1.1 mjy for intrinsic pulse 
widths of 10°-90°. 

The signal-to-noise ratio used in the sensitivity equation 
is calculated in the time domain. We note that some pulsar 
surveys (e.g. Camilo, Nice & Taylor 1996; Manchester et al. 
2001) have based their sensitivity analysis on the so-called 
'spectral' S/N, computed from the amplitudes of harmonics 
in the power spectrum. Our search code computes such a 
value for all candidates, however this is used only for the se- 
lection of signals to be subjected to a fine search in the time 
domain, and the threshold value used (5.0) is sufficiently 
low that we expect no loss of significant candidates at this 
stage. The (time domain) S/N threshold of 8.0 imposed in 
the selection of candidates for human viewing is also irrel- 
evant for the sensitivity analysis since it found that nearly 
all candidates re-detected in subsequent observations had 
htmS/Ns greater than 10.0, despite attempts to re-detect large 
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Figure 7. Histograms depicting distribution in pulse period of 
new pulsars (dashed line) and previously known pulsars in the 
survey region (solid line). 



numbers of candidates with S/N between 8.0 and 10.0. We 
therefore use a = 10.0 in the above analysis and in the 
curves plotted in Figure ^[ 

We caution that this analysis should only be taken as 
approximate. The variation of a pulsar's flux due to scintil- 
lation adds a considerable degree of uncertainty to its de- 
tectability, and particularly for millisecond pulsars lack of 
time resolution and prevalence of interference signals cause 
the human viewer to effectively adopt a higher threshold 
S/N. In particular, whilst this analysis suggests that this 
and similar surveys could be significantly sensitive to sub- 
millisecond pulsars, we would treat such a claim with skep- 
ticism. A sub-millisecond pulsar with P = 0.8 ms and a 
moderately low dispersion measure of 25 cm -3 pc would ex- 
perience 0.25 ms of dispersion smearing, resulting in a pulse 
profile with a width of at least 110°. Combined with the 
fact that only 4 or perhaps 8 bins would be present in the 
pulse profile, the human viewer is forced to judge the ori- 
gin of the signal essentially entirely upon the reported S/N, 
which itself becomes subject to significant uncertainty when 
the number of bins is so few. We therefore expect that the 
standard sensitivity analysis underestimates the minimum 
detectable flux density by a factor of a few for pulsars with 
periods shorter than a millisecond. Likewise, the minimum 
detectable flux density may be underestimated for very slow 
(P > 5 s) pulsars due to the rather short time constant 
(0.9 s) employed in the digitiser system. The P -1 / 2 duty 
cycle dependence of slow pulsars (e.g. Rankin 1990) aids 
the situation somewhat, but nevertheless it is expected that 
a few percent of very slow pulsars have pulse widths greater 
than 0.5 s, and these would experience significant S/N loss. 

A total of 101 previously known pulsars were re- 
detected in the survey. Due to the large angular extent of 
each of the tessellated units of four pointings, the actual sur- 
vey regions are significantly non-rectangular in / and b (see 
Figure |2|) and three of these pulsars actually lie outside the 
nominal survey region. There were 135 previously discov- 




50 100 
Signal to noise ratio 

Figure 8. Histograms depicting distribution in maximum de- 
tected signal to noise ratio (S/N) of new pulsars (dashed line) 
and previously known pulsars in the survey region (solid line). 
One new and 8 known pulsars had S/Ns greater than 120 and are 
not included in this plot for clarity in low S/N bins. 
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Figure 9. Histograms depicting distribution in dispersion mea- 
sure of new pulsars (dashed line) and previously known pulsars 
in the survey region (solid line). 



ered pulsars in the declared region, leaving 37 undetected in 
the survey. Tables jj] and |^ list the detected and undetected 
known pulsars, including for each pulsar the period, dis- 
persion measure, flux density at 1400 MHz (where known), 
position in Galactic coordinates, angular offset from the cen- 
tre of the beam in which it was detected (where known, or 
the nearest processed beam for undetected pulsars) , and the 
signal to noise ratio of detections. In most cases the values 
for period, DM and flux density derive from the works of 
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Table 3. Dispersion measure, astrometric and spin parameters for new slow pulsars 

Name a (J2000) S (J2000) P P Epoch P DM 

(s) (MJD) (lCT 15 ) (cm" 3 pc) 



J0834- 


-60 


08 h 34 m 50(40) 


-60°35(5)' 


0.384645(6) 


51401.1 




20(6) 


J0843- 


-5022 


08 h 43 m 09?884(8) 


-50°22'43'.'10(8) 


0.2089556931527(14) 


51500.0 


0.17238(14) 


178.47(9) 


J0849- 


-6322 


08 h 49 m 42f59(2) 


-63°22'35'.'0(1) 


0.367953256307(5) 


51500.0 


0.7908(5) 


91.29(9) 


J0932- 


-3217 


09 h 32 m 39?15(6) 


-32°17'14'.'2(8) 


1.93162674308(19) 


51500.0 


0.250(15) 


102.1(8) 


J0934- 


-4154 


09 h 34 m 58?20(3) 


-41°54'19'.'5(3) 


0.570409236430(14) 


51650.0 


0.269(3) 


113.79(16) 


J1055- 


-6905 


10 h 55 m 44!71(9) 


-69° 05' 11'.' 4 (4) 


2.9193969868(3) 


51500.0 


20.336(15) 


142.8(4) 


J1057- 


-47 


10 h 57 m 45(30) 


-47° 57(5)' 


0.62830(3) 


50989.1 




60(8) 


J1204- 


-6843 


12 h 04 m 36!72(l) 


-68°43'17'.'19(8) 


0.3088608620097(19) 


51500.0 


0.21708(19) 


133.93(9) 


J1215- 


-5328 


12 h 15 m 0O?62(7) 


-53°28'31'.'6(7) 


0.63641413680(5) 


51500.0 


0.115(4) 


163.0(5) 


J1231- 


-47 


12 h 31 m 40(30) 


-47° 46(5)' 


1.8732(3) 


51200.8 




31(30) 


J1236- 


-5033 


12 h 36 m 59?15(l) 


-50°33'36'.'3(1) 


0.294759771191(4) 


51500.0 


0.1556(4) 


105.02(11) 


J1244- 


-5053 


12 h 44 m ll?48(l) 


-50°53'20'.'6(1) 


0.275207111323(4) 


51500.0 


0.9998(4) 


109.95(12) 


J1352- 


-6803 


13 h 52 m 34?45(4) 


-68°03'37'.'1(4) 


0.628902546380(16) 


51650.0 


1.234(3) 


214.6(2) 


J1410- 


-7404 


14 h 10 m 07?370(5) 


-74°04'53'.'32(2) 


0.2787294436271(15) 


51460.0 


0.00674(9) 


54.24(6) 


J1414- 


-6802 


14 h 14 m 25!7(l) 


-68°02'58(1)" 


4.6301880619(4) 


51650.0 


6.39(7) 


153.5(6) 


J1415- 


66 


14 h 15 m 25(50) 


-66° 19(5)' 


0.392480(10) 


51396.2 




261(6) 


J1424- 


-69 


14 h 24 m 15(60) 


-69° 56(5)' 


0.333415(8) 


51309.7 




123(4) 


J1517- 


-4356 


15 h 17 m 27!34(l) 


-43°56'17'.'9(2) 


0.650836871901(6) 


51500.0 


0.2155(6) 


87.78(12) 


J1528- 


-4109 


15 h 28 m 08?033(8) 


-41°09'28'.'8(2) 


0.526556139140(4) 


51500.0 


0.3955(4) 


89.50(10) 


J1531- 


-4012 


15 h 31 m 08!05(l) 


-40°12'30'.'9(4) 


0.356849312855(6) 


51500.0 


0.0963(6) 


106.65(12) 


J1535- 


-4114 


15 h 35 m 17!07(l) 


-41°14'03'.'1(3) 


0.432866133845(6) 


51500.0 


4.0705(6) 


66.28(14) 


J1536- 


-44 


15 h 36 m 15(30) 


-44° 16(5)' 


0.46842(6) 


51063.2 




110(30) 


J1537- 


49 


15 h 37 m 30(30) 


-49°09(5)' 


0.301313(6) 


51402.3 




65(4) 


J1540- 


-63 


15 h 40 m 20(40) 


-63° 24(5)' 


1.63080(16) 


51307.7 




160(20) 


J1603- 


-3539 


16 h 03 m 53?697(5) 


-35°39'57'.'1(3) 


0.1419085889640(9) 


51650.0 


0.12425(17) 


77.5(4) 


J1617- 


-4216 


16 h 17 m 23?38(5) 


-42°16'59(1)" 


3.42846630955(13) 


51500.0 


18.129(15) 


163.6(5) 


J1641- 


-2347 


16 h 41 m 18!04(6) 


-23°47'36(6)" 


1.091008429855(16) 


51500.0 


0.0411(15) 


27.7(3) 


J1649- 


-5553 


16 h 49 m 31!l(l) 


-55°53'40(2)" 


0.61357070436(7) 


51650.0 


1.698(16) 


180.4(11) 


J1655- 


-3048 


16 h 55 m 24!53(2) 


-30°48'42(1)" 


0.542935874228(9) 


51500.0 


0.0366(9) 


154.3(3) 


J1701- 


-3130 


17 h 01 m 43!513(5) 


-31°30'36'.'7(4) 


0.2913414710251(12) 


51500.0 


0.05596(12) 


130.73(6) 



Taylor, Manchester & Lyne (1993), Lorimer et al. (1995) or 
D'Amico et al. (1998). 

From the tabulated position offsets to the nearest pro- 
cessed beams and given the fact that the centres of adjacent 
beams of the grid are spaced 14 arcminutes apart on the sky, 
it is clear that the reason for many of the non-detections 
was incomplete sky coverage. Twelve undetected pulsars lay 
greater than 10 arcminutes from the nearest beam, leading 
to an alternative completeness metric of 1 — 12/135 = 91 
per cent, (or 1 — 9/135 = 93 per cent if the loss is off- 
set by the three known pulsars detected outside the survey 
region). Of the remaining 25 non-detections, 11 have pub- 
lished flux densities near 1400 MHz and are plotted along 
with the sensitivity curves of Figure ^. Allowing for scin- 
tillation, all are compatible with having flux densities be- 
low the sensitivity limit. We therefore expect that most of 
those pulsars lacking flux density measurements were also 
below the detection threshold at the time of observation, 
and conclude that the search procedure was adequate and 
robust in its rejection of interference without significant loss 
of pulsars. Possible exceptions to this statement are PSR 
B1556-44 which was detected at 17 times its true spin fre- 
quency (probably due to proximity in period to a persistent 



256- ms interference signal), and perhaps PSR J 1130-6807 
which has a similar pulse period but being of unknown flux 
density may have simply fallen below the sensitivity limit. 
The discovery of PSR J1712-2715 (below) with P = 255.4 
ms indicates that rough proximity to interference signals 
is not always problematic. The millisecond pulsar (MSP) 
fraction of the undetected pulsars is high, however all ex- 
cept J1911-1114 (Lorimer et al. 1996) were discovered in 
deep directed searches of globular clusters (Lyne et al. 1987; 
Biggs et al. 1994; D'Amico et al. 2001) which found mainly 
millisecond pulsars. Figure ^] shows a histogram of flux den- 
sities for previously known pulsars of published flux density 
with processed beams centred less than 10' away, with the 
distribution of non-detections hashed. As expected from the 
sensitivity curves shown in Figure ^, the survey was sensitive 
to most pulsars brighter than 1 mjy, insensitive to pulsars 
with S < 0.1 mjy, and recorded a mixture of detections and 
non-detections in the remaining range due to scintillation 
and the distribution of pulse widths. 

Examination of the detection parameters of previously 
known pulsars reveals some discrepancies with previously 
published results. From inspection of the position offsets of 
newly discovered pulsars (below, Table ^) , we estimate a po- 
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Table 3 - continued 

Name a (J2000) <5 (J2000) P P Epoch P DM 

(s) (MJD) (lCT 15 ) (cm- 3 pc) 



J1706-61 
J1712-2715 
J 171 7-5800 
J 172 1-1939 
J1739-1313 

J1741-2019 
J1742-4616 
J1743-4212 
J1744-1610 
J1745-0129 

J1802+0128 

J1805-0619 

J1806+10 

J1808-3249 

J1809-0743 

J1811-0154 

J1819+1305 

J1824-25 

J1832-28 

J1837+1221 

J1837-1837 

J1842+1332 

J1848+12 

J1855-0941 

J1857-1027 

J1901-1740 
J1919+0134 
J1943+0609 
J1947+0915 
J1956+0838 



17 h 06 m 40(40) 

17 h 12 m ll?71(l) 

17 h 17 m 35?65(2) 

17 h 21 m 46?61(4) 

17 h 39 m 57!821(6) 

17 h 41 m 06?87(3) 

17 h 42 m 26?10(2) 

17 h 43 m 05! 223(5) 
17 h 44 m 16 s 5 3 4 ( 7 ) 

17 h 45 m 02f06(l) 

18 h 02 m 27!45(2) 

18 h 05 m 31?436(9) 

18 h 06 m 50(20) 

18 h 08 m 04f48(2) 

18 h 09 m 35f92(l) 

18 h ll m 19!88(3) 
18 h 19 m 56?22(4) 
18 h 24 m 15(20) 
18 h 32 m 30(20) 
18 h 37 m 07fl2(4) 

18 h 37 m 54?25(l) 

18 h 42 m 29!96(6) 

18 h 48 m 30(20) 

18 h 55 m 15!68(3) 

18 h 57 m 26?45(5) 

19 h 01 m 18!03(6) 
19 h 19 m 43?62(3) 
19 h 43 m 29! 132(5) 
19 h 47 m 46?22(5) 
19 h 56 m 52?26(2) 



-61°11(5)' 

-27°15'53(2)" 

-58°00'05"4(3) 

-19°39'49(5)" 

-13°13'18('6(4) 

-20°19'24(5)" 

-46°16'53('5(4) 

-42°12'02^'4(2) 

-16°10'35^'8(8) 

-01°29'18^'1(4) 

+01°28'23'.'7(4) 

-06°19'45^'4(4) 

+10°24(5)' 

-32°49'34(1)" 

-07°43'01"4(5) 

-01°54'30'.'9(7) 
+13°05'14'.'2(7) 
-25°36(5)' 
-28° 43(5)' 
+12°21'54'.'0(6) 

-18°37'08(2)" 

+13°32'01'.'5(9) 

+12°50(5)' 

-09°41'02(1)" 

-10°27'01(2)" 

-17°40'00(6)" 

+01°34'56'.'5(7) 

+06°09'57'.'6(1) 

+09°15'08'.'0(8) 

+08°38'16'.'8(4) 



J2007+0809 20 h 07 m 13?5(l) +08°09'33(2)" 0.32572436605(5) 



0.361922(8) 51308.7 

0.255359660118(3) 51500.0 

0.321793346869(6) 51650.0 

0.404039751280(15) 51500.0 

1.215697613611(9) 51500.0 

3.90450636119(13) 51500.0 

0.412401047219(7) 51650.0 

0.3061669878595(16) 51650.0 

1.757205868816(16) 51500.0 

1.045406855598(18) 51650.0 

0.554261603931(10) 51650.0 

0.454650713078(7) 51650.0 

0.484285(15) 51259.8 

0.364912241765(10) 51500.0 

0.313885674748(5) 51650.0 

0.92494482303(4) 51650.0 

1.06036354400(6) 51650.0 

0.223319(3) 51067.5 

0.199300(3) 51064.3 

1.96353198352(12) 51650.0 

0.618357697387(16) 51500.0 

0.47160357893(3) 51650.0 

0.75473(7) 51316.7 

0.34540115992(4) 51500.0 

3.6872190477(3) 51650.0 

1.95685759005(16) 51500.0 

1.60398355528(6) 51650.0 

0.446226281658(3) 51650.0 

1.48074382424(9) 51650.0 

0.303910924347(7) 51650.0 

51650.0 





78(6) 


1.2793(3) 


92.64(13) 


0.1957(10) 


125.22(14) 


0.1283(15) 


103(2) 


0.0817(9) 


58.2(5) 


16.260(13) 


74.9(4) 


0.0338(12) 


115.96(14) 


0.7834(4) 


131.94(5) 


2.3767(16) 


66.67(14) 


0.631(4) 


90.1(11) 


2.109(3) 


97.97(12) 


0.9690(13) 


146.22(9) 




58(6) 


7.0494(10) 


147.37(19) 


0.1521(9) 


240.70(14) 


1.608(6) 


148.1(3) 


0.373(9) 


64.9(4) 




155(3) 




127(3) 


6.200(16) 


100.6(4) 


5.4950(12) 


100.74(13) 


0.229(7) 


102.5(7) 




139(20) 


0.240(3) 


152.2(3) 


10.55(6) 


108.9(7) 


0.823(16) 


24.4(6) 


0.589(11) 


191.9(4) 


0.4659(6) 


70.76(6) 


0.478(16) 


94(4) 


0.2199(13) 


68.2(13) 


0.137(7) 


53.9(10) 



sition uncertainty of 7' for detections in this survey, consis- 
tent with the beam spacing. Of those detected pulsars with 
published positions greater than 10' from the beam centre, 
two were not detected in pointings made closer to the pub- 
lished position. The (B1950) declinations reported for PSRs 
B1232-55 and B1359-51 of -55°00(10)' and -51° 10(15)' 
(Manchester et al. 1978) respectively are inconsistent with 
the B1950 declinations of -54° 40(5)' and -50°06(5)' of the 
detections made in this survey. The published right ascen- 
sion values have much smaller uncertainties (due to the 
shape of the beam of the Molonglo telescope with which 
they were discovered), and are consistent with our detec- 
tions. 

As was spectacularly illustrated with PSR J2144-3933 
(Young, Manchester & Johnston 1999), pulsar surveys in the 
past have been prone to detecting pulsars with an apparent 
pulse frequency an integer multiple of the true frequency. 
We expect that this is due to the flattening of the power 
spectrum with a boxcar filter, the presence of interference 
signals of similar period, the dominance of odd harmonics in 
the pulse shape, or the exclusion of the fundamental as being 
below a cutoff frequency (as was the case for J2144-3933). 



We found that PSRs J1403-7646, J1817-3837 and J1901- 
0906 (Lyne et al. 1998) actually possess periods a factor 
of two greater than the published values. Conversely, pul- 
sars J1036-4926, B1110-69, B1524-39, B1556-44, B1706- 
16, B1709-15, B1717-16, B1718-19 and B1848+12 were er- 
roneously re-detected with shorter pulse periods. All but 
one of these results were made in the early stages of pro- 
cessing before interference rejection had been fine-tuned to 
avoid accidental filtering of pulsar harmonics near interfer- 
ence signals. After correction for such errors and with the 
exception of the three pulsars listed above, all other detec- 
tions were made with pulse periods consistent with previ- 
ously published parameters. 

3.2 Newly Discovered Pulsars 

The initial viewing of the survey results and subsequent con- 
firmation observations resulted in the discovery of 58 new 
pulsars, 8 of which possess short spin periods and small pe- 
riod derivatives indicative of recycling and have been re- 
ported elsewhere (Edwards & Bailes 2001b; Edwards & 
Bailes 2001a). The final careful review of candidates pro- 



© 0000 RAS, MNRAS 000, 000-000 



14 R. T. Edwards, M. Bailes, W. van Straten and M. C. Britton 



Table 4. Detection parameters, pulse widths and derived parameters for new slow pulsars 

Name A pos S/N w$q wio I b d \z\ t c B E 

(') (°) (°) (°) (°) (kpc) (kpc) (Myr) (10 12 G) (10 30 erg s" 1 ) 



J0834- 


-60 




13.9 


J0843- 


-5022 




10.1 


J0849- 


-6322 


3.9 


13.2 


J0932- 


-3217 


5.9 


17.1 


J0934- 


-4154 


4.4 


10.4 


J1055- 


-6905 


3.5 


17.8 


J1057- 


-47 




16.3 


J1204- 


-6843 


5.5 


18.5 


J1215- 


-5328 


8.1 


12.9 


J1231- 


-47 




28.5 


J1236- 


-5033 


3.2 


13.2 


J1244- 


-5053 


4.3 


12.4 


J1352- 


-6803 


3.5 


32.7 


J1410- 


-7404 


5.7 


30.4 


J1414- 


-6802 


2.0 


27.3 


J1415- 


66 




26.3 


J1424- 


-69 




13.0 


J1517- 


-4356 


5.6 


11.8 


J1528- 


-4109 


3.7 


19.8 


J1531- 


-4012 


5.3 


15.4 


J1535- 


-4114 


3.8 


95.0 


J1536- 


-44 




15.8 


J1537- 


19 




14.5 


J1540- 


-63 




25.3 


J1603- 


-3539 


4.1 


13.3 


J1617- 


-4216 


5.1 


10.3 


J1641- 


-2347 


3.7 


85.1 


J1649- 


-5553 


7.3 


13.5 


J1655- 


-3048 




32.3 


J1701- 


-3130 


5.9 


26.8 



-83.9 



11.0 


50.0 


-91.5 


7.0 


166.3 


-80.6 


3.4 


7.6 


-98.7 


10.1 




-91.6 


5.8 


10.7 


-67.1 


11.5 




-76.0 


10.6 


15.2 


-61.3 


22.6 




-62.5 


30.7 




-60.5 


14.7 


20.8 


-59.4 


8.1 




-58.2 


16.0 


188.9 


-51.4 


2.3 


4.5 


-51.7 


7.8 


14.0 


-49.4 


6.5 




-48.8 


3.6 




-49.2 


6.1 


14.4 


-31.1 


6.4 


12.3 


-27.9 


11.1 




-26.9 


12.5 


18.3 


-26.8 






-28.5 


26.3 




-31.3 


12.4 




-39.4 


30.4 




-18.8 


3.7 


9.5 


-21.5 


16.2 


29.7 


-4.2 


39.9 




-27.9 


57.0 


72.9 


-7.8 


11.6 


23.5 


-7.5 



-11.9 


0.5 


0.10 


-4.9 


7.7 


0.66 


-12.2 


> 8.4 


> 1.8 


14.1 


3.8 


0.93 


7.4 


3.2 


0.41 



-8.5 


> 12 


> 1.8 


10.7 


3.0 


0.56 


-6.2 


5.7 


0.61 


9.0 


> 11 


> 1.8 


15.0 


1.6 


0.42 


12.2 


> 8.3 


> 1.8 


12.0 


> 8.5 


> 1.8 


-5.9 


14 


1.4 


-12.0 


2.1 


0.45 


-6.4 


6.6 


0.74 


-4.8 


14 


1.2 


-8.5 


6.3 


0.94 


11.5 


4.4 


0.88 


12.7 


6.0 


1.3 


13.1 


> 7.8 


> 1.8 


11.8 


2.8 


0.57 


9.3 


5.3 


0.85 


5.2 


1.7 


0.16 


-6.5 


7.5 


0.84 


12.5 


3.8 


0.83 


5.9 


6.3 


0.65 


14.7 


1.3 


0.34 


-7.1 


14 


1.7 


7.9 


8.6 


1.2 


6.1 


4.5 


0.50 



19.2 


0.19 


746 


7.37 


0.55 


627 


122 


0.70 


1.37 


33.5 


0.40 


57.3 


2.27 


7.8 


32.3 


22.5 


0.26 


291 


87.4 


0.27 


17.7 


30.0 


0.22 


240 


4.36 


0.53 


1894 


8.07 


0.89 


196 


655 


0.044 


12.3 


11.5 


5.5 


2.54 


47.8 


0.38 


30.9 


21.1 


0.46 


107 


58.7 


0.19 


83.6 


1.68 


1.3 


1981 


18.1 


0.13 


1716 


3.00 


8.0 


17.8 


421 


0.21 


1.25 


5.73 


1.0 


290 


235 


0.14 


9.02 


82.5 


0.13 


89.3 



duced a further 11 slow pulsars. Data from all pulsars were 
folded at twice and three times the discovery period to de- 
tect any errors of the kind described in the previous section. 
Pulsars J1055-6905, J1517-4356, J1802+0128 and J1808- 
3249 were all initially discovered at half the true spin period, 
however J1517-4356 and J1808-3249 were later detected at 
the correct period in subsequent survey observations, with 
higher signal to noise ratio. Pulse profiles for the 61 slow 
pulsars discovered in the survey are presented in Figure ^| 
For those pulsars with a timing solution the profile arises 
from the summation of numerous good observations, whilst 
for others the profile from the single best observation made 
to date is provided. It should be noted that the baselines 
of some profiles are corrupted due to the response of the 
filterbank/digitiser low pass filter. 

Timing solutions have been obtained for 48 slow pul- 
sars. The data for each pulsar span 490-830 days, depend- 
ing on when the candidate re-observation was made. Since 
all pulsars have been timed for well over a year, we were 
able to accurately measure positions and period derivatives, 
in addition to making improved measurements of periods 
and dispersion measures. We have not seen obvious timing 
noise in any pulsar, however it is possible that such noise 



is present and is absorbed by deviations from the "true" 
values of the fitted parameters. Dispersion measures were 
fitted for by the use of TOAs derived from sub-divisions of 
the observing band, or where available, 660 MHz observa- 
tions. In the latter case the signal to noise ratio available 
in a reasonable integration time was poor and the use of an 
independent template profile was eschewed in favour of the 
superior 20-cm template. 

Table ^ presents the basic parameters of the slow pul- 
sars discovered in this survey, derived from timing measure- 
ments when available. Values in parentheses denote uncer- 
tainties in the last quoted digit and represent twice the for- 
mal uncertainty produced by TEMPO. In the case of pulsars 
lacking a timing solution we use twice the errors derived 
from the fine P-DM search, or for positions, errors corre- 
sponding to a randomly oriented offset of 7 arcminutes (i.e. 
7'/V2 ~ 5' in each of a and 5). PSR J1802+0128 was timed 
for 493 days before it was discovered that the true spin pe- 
riod was twice the assumed value. Values presented for P 
and P are twice the values derived by this timing analy- 
sis and retain their prior relative errors. Due to the effects 
of the evolution of profile morphology with frequency and 
the use of the 20 cm as template for timing 50 cm observa- 
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Table 4 - continued 

Name A pos S/N w$o W W I b d \z\ t c B E 

(') (°) (°) (°) (°) (kpc) (kpc) (Myr) (10 12 G) (10 30 erg s" 1 ) 



J1706-61 




21.6 


2.8 




-30.8 


-12.1 


3.8 


0.80 








J1712-2715 


1.1 


41.7 


46.5 


90.2 


-2.7 


7.1 


3.1 


0.38 


3.16 


0.58 


3033 


J 171 7-5800 


5.0 


10.9 


23.5 




-27.3 


-11.5 


> 8.8 


> 1.8 


26.1 


0.25 


232 


J1721-1939 


4.6 


9.7 


16.9 




4.9 


9.6 


4.7 


0.78 


49.9 


0.23 


76.8 


J1739-1313 


3.3 


22.6 


2.4 


5.3 


12.8 


9.3 


2.0 


0.33 


236 


0.32 


1.80 


J1741-2019 


6.0 


59.6 


11.1 


13.8 


6.8 


5.4 


2.0 


0.19 


3.80 


8.1 


10.8 


J1742-4616 


3.1 


26.3 


21.8 


33.2 


-15.2 


-8.5 


5.0 


0.74 


193 


0.12 


19.0 


J1743-4212 


5.2 


19.6 


8.6 


18.8 


-11.6 


-6.5 


4.7 


0.53 


6.19 


0.50 


1078 


J1744-1610 




12.8 


7.4 


10.2 


10.8 


6.9 


2.0 


0.24 


11.7 


2.1 


17.3 


J1745-0129 


5.4 


12.2 


2.4 


8.3 


23.8 


14.0 


> 7.3 


> 1.8 


26.3 


0.82 


21.8 


J1802+0128 


1.9 


12.5 


1.1 




28.6 


11.6 


> 8.8 


> 1.8 


4.16 


1.1 


489 


J1805-0619 


6.0 


11.1 


12.7 


23.0 


22.0 


7.2 


6.7 


0.84 


7.43 


0.67 


107 


J1806+10 




66.3 


6.7 




37.3 


14.6 


1.0 


1.0 








J1808-3249 


2.1 


40.8 


13.8 


20.7 


-1.0 


-6.1 


5.1 


0.54 


0.82 


1.6 


5727 


J1809-0743 


1.3 


13.0 


13.8 




21.2 


5.7 


> 18 


> 1.8 


32.7 


0.22 


194 


J1811-0154 


3.8 


21.8 


8.2 


77.8 


26.6 


8.0 


9.6 


1.3 


9.12 


1.2 


80.2 


J1819+1305 


5.1 


14.3 


21.4 




41.2 


12.8 


4.4 


0.98 


45.0 


0.64 


12.4 


J1824-25 




15.6 


9.3 




7.1 


-5.9 


5.2 


0.53 








J1832-28 




18.0 


21.3 




5.1 


-8.9 


6.4 


1.00 








J1837+1221 


4.9 


16.1 


3.1 




42.4 


8.7 


6.1 


0.93 


5.02 


3.5 


32.3 


J1837-1837 


4.8 


11.0 


5.8 


11.6 


14.8 


-5.5 


3.0 


0.29 


1.78 


1.9 


918 


J1842+1332 


3.5 


33.1 


72.5 




44.1 


8.1 


5.9 


0.83 


32.6 


0.33 


86.2 


J1848+12 




17.8 


51.0 




44.1 


6.5 


7.5 


0.84 








J1855-0941 


2.9 


19.1 


26.1 




24.7 


-5.2 


4.9 


0.45 


22.8 


0.29 
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J1857-1027 


0.9 


71.4 


14.6 


20.5 


24.3 


-6.1 


3.6 


0.38 


5.54 


6.3 


8.31 


J1901-1740 


5.5 


131.0 


6.6 


22.1 


18.1 


-10.1 


1.3 


0.22 


37.7 


1.3 


4.34 


J1919+0134 


5.0 


34.2 


10.5 


18.2 


37.6 


-5.6 


10 


0.99 


43.1 


0.98 


5.64 


J1943+0609 


2.5 


30.1 


8.9 


17.7 


44.5 


-8.6 


3.9 


0.58 


15.2 


0.46 
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J1947+0915 


6.1 


12.7 


7.5 




47.7 


-8.1 


5.8 


0.81 


49.1 


0.85 


5.81 


J1956+0838 


5.8 


15.0 


15.4 




48.3 


-10.3 


1.3 


0.77 


21.9 


0.26 


309 


J2007+0809 


5.2 


14.4 


92.8 




49.2 


-12.8 


3.4 


0.76 


37.7 


0.21 


156 



tions, the formal errors presented for the dispersion measure 
may have been underestimated. Pulsar names are assigned 
from their equatorial coordinates in the J2000.0 equinox, 
with four digits of declination for those with accurate posi- 
tions from timing solutions, and two digits for those with- 
out. The latter names are provisional and will be altered 
when accurate timing positions become available. Table ^ 
lists additional parameters for the newly discovered pulsars, 
including the best S/N of the discovery observation (s) and 
the pulsar's position offset from the centre of the beam, the 
width of the pulse profile at 10 and 50 per cent of peak in- 
tensity, its position in Galactic coordinates, its distance and 
^-displacement from the Galactic plane under the model of 
Taylor & Cordes (1993) (accurate to 30 per cent on aver- 
age), and inferred parameters concerning the pulsar spin- 
down. These assume magnetic dipole spin-down and com- 
prise the characteristic age (t c = P/2P), surface magnetic 

field strength (B = 3.2 x 10 19 G s~ 1/2 \/pP) and spin-down 
power (E — 4ir 2 IPP~ S , assuming I — 10 45 g cm 2 for the 
moment of inertia of the neutron star). Both tables are ac- 
cessible on the internet in machine-readable format at the 



Swinburne Pulsar Group home pageQ For several pulsars the 
dispersion measure is higher than that allowed in the given 
direction under the model of Taylor & Cordes (1993), and 
the values presented are given as lower limits. Such findings 
are not uncommon (e.g. Camilo & Nice 1995; D'Amico et al. 
1998), and indicate that the model probably underestimates 
the scale height of the Galactic electron distribution. 

The spin parameters of the newly discovered systems 
are similar to those of pulsars previously known in the search 
region. Figure [] shows the distribution in period and period 
derivative of new pulsars with timing solutions and of pre- 
viously known pulsars inside and outside the survey region. 
Both new and previously known slow pulsars in the region 
tend to have longer inferred characteristic ages than those 
outside the survey region, by simple virtue of the fact that 
pulsars are born near the Galactic plane and typically take 
several Myr to reach a z-height corresponding to |fe| > 5° (for 
typical distances of several kpc). Figure ^ shows histograms 
of pulse period for the new and previously known popu- 

t http:/ /www. astronomy.swin.edu.au/pulsar 
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Figure 10. Histograms depicting distribution in angular dis- 
placement from the Galactic plane of new pulsars (dashed line) 
and previously known pulsars in the survey region (solid line). 



lation. It appears that this survey has uncovered a higher 
fraction of pulsars in the period range of 6-50 ms, however 
this effect is not highly significant : a Kolmogorov-Smirnov 
(K-S) test on the distributions for P < 100 ms yields a 
46 per cent probability of the two samples arising from the 
same parent distribution. When binary evolution considera- 
tions are taken into account and the new sample restricted to 
the four MSPs with probable helium companions the signif- 
icance rises to 91 per cent, however due to the small sample 
size this result must be treated with caution (Edwards & 
Bailes 2001a). 

The discrepancy in the case of the slow pulsar popu- 
lation is more significant, with a K-S test indicating a dif- 
ferent distribution at the 98 per cent level. The reason for 
the discrepancy, which is mainly seen as a deficit of pulsars 
with P ~ 1 s, is not well understood. The problem remains 
(at 94 per cent significance) when the new sample is com- 
pared only to those previously known pulsars re-detected in 
this survey, confirming that the survey was able to detect 
pulsars in the period range in question. However, if there 
existed (by an unknown mechanism) a reduced sensitivity 
(or higher effective threshold signal to noise ratio) around 
P ~ 1 s, one might expect the period distribution of the new 
pulsars to be more strongly affected than that of the previ- 
ously known pulsars since the new population is on average 
of lower flux density. We note that the rejection of pulsars 
as mis-categorized interference signals cannot explain this 
result, since in this case one would expect an equal rate of 
rejections of new and previously known pulsars independent 
of signal to noise ratio or flux density. Comparison of the 
(yet to be measured) flux densities of new pulsars in and 
out of the depleted period range will help in evaluating the 
effective sensitivity of the survey as a function of period. 

Figure ^ shows the distribution of signal to noise ratios 
in the best detections of new and previously known pulsars. 
As one would expect, most of the pulsars with high signal 
to noise ratios were detected in earlier surveys. A few of the 



new pulsars however were very strong and may have been 
missed in earlier surveys due to scintillation or due to an in- 
trinsically flat spectrum. It is apparent from the histogram 
and from inspection of Tables l] and |^ that the threshold 
signal to noise ratio for this survey is approximately 10, in 
contrast to the value of 8.0 commonly used in previous sur- 
veys in assessing sensitivity. Numerous promising candidates 
with signal to noise ratios in the range of 8-10 were sub- 
jected to re-observation however only one was re-detected 
in such observations despite using longer integration times, 
probably because they actually arose from interference or 
by random chance (given the size and dimensionality of the 
search space). 

As noted in the introduction, high frequency surveys are 
also expected to sample a different area of the pulsar spectral 
index distribution, compared to low frequency surveys. The 
sensitivity of the present work to a pulsar with a spectral 
index of —1.7 (typical of those discovered at 70 cm; Toscano 
et al. 1998) is comparable to that of the Parkes Southern 
Pulsar Survey (Manchester et al. 1996), the most sensitive 
previous 70 cm search to cover a large region of the area ob- 
served by the present study. One therefore expects the dis- 
tances of the newly discovered pulsars to be comparable to 
the previously known population in the region, and as shown 
in Figure^, this is indeed the case. The bulk of newly discov- 
ered pulsars had S/Ns ^ 30, suggesting (in conjunction with 
their non-detection in the 70 cm survey) spectral indices of 
up to —0.7. Eleven newly discovered pulsars are visible from 
the 305-m Arecibo telescope and were presumably within 
the search area of previous surveys conducted there (Foster 
et al. 1995; Camilo et al. 1996; Ray et al. 1996; Lommen 
et al. 2000), which were typically sensitive to (S — 0.1) pul- 
sars brighter than ~ 1 mjy at 70 cm. Their non-detection in 
the Arecibo surveys suggests a J> —1.1. Several pulsars were 
discovered at high S/N, which could be indicative of positive 
spectral indices (for example, J1806+10 in Arecibo territory, 
yielding a ~ 2.6), however the non-detections may well be 
the result of incomplete surveys or interstellar scintillation. 

Accurate characterization of the spectral indices of 
the newly discovered pulsars must await calibrated multi- 
frequency flux density measurements, however we note in 
passing that the detection of these pulsars at 660 MHz re- 
quired significantly more integration time than expected for 
pulsars of average spectral index. As indicated by the lack 
of any enhanced preference for low Galactic latitudes in the 
newly detected sample compared to the previously known 
pulsars (Figure [lo|), reduced scatter-broadening at 20-cm in 
general does not appear to have been a significant factor in 
the discovery of the new pulsars. 

3.3 Individual Pulsars of Interest 

As previously reported (Edwards & Bailes 2001a), the pulse 
profile of PSR J1410-7404 is exceedingly narrow and con- 
tradicts the pulse width - period relation of Rankin (1990). 
Since all major contradictions in the past have been from ap- 
parently recycled pulsars, it is conceivable that J1410-7404 
is also recycled, a hypothesis supported by the small mag- 
netic field strength inferred from timing observations. The 
newly discovered pulsar J 1706-61 also has a measured pro- 
file width seemingly in disagreement with Rankin (1990), in 
this case being 70 per cent of the predicted minimum width. 
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Figure 11. Greyscale plots of detected flux density as a function of pulse phase and pulse number for four newly discovered pulsars. 
Each row represents a single pulse. 



However, we caution that this measurement derives from a 
single observation of moderate signal-to-noise ratio (as visi- 
ble in Figure and that a more definite conclusion awaits 
the availability of an extended data set from ongoing tim- 
ing observations. Should the improved profile maintain the 
narrow width derived here, the magnetic field strength and 
characteristic age derived from timing measurements will 
be of great utility in evaluating the recycling hypothesis for 
PSRs J1410-7404 and J1706-61. 

As expected from the relatively old age of the sample 
detected here, numerous pulsars appear to exhibit noticeable 
pulse nulling (e.g. Ritchings 1976) . Several pulsars also show 
drifting sub-pulses (e.g. Backer 1973) of varying degrees of 
regularity. Figure [ill depicts the instantaneous flux density 
of four newly discovered pulsars as a function of pulse phase 
and pulse number. It is apparent that the emission of PSRs 
J1231-47 and J1919+0134 is strongly modulated by sub- 
pulses showing very regular drift. PSR J1857-1027 also ap- 
pears to show drifting sub-pulses in this and other observa- 
tions, although as a result of the short duration of observa- 
tions made for timing analysis, combined with the long pulse 
period (~ 3.7 s), the number of pulses recorded and hence 
the conspicuousness of this effect in archival observations is 



reduced. Detailed analysis of these pulsars will appear in a 
forthcoming paper (Ord et al. in preparation). 

Pulse nulling appears to occur on a wide range of time 
scales in the detected pulsars. The data for PSR J1231-47 
in Figure [ll] were recorded in a second attempt at confir- 
mation of this pulsar, resulting in a signal to noise ratio of 
82 in an observation of 530 s. Several subsequent observa- 
tions of between 900 and 4000 s succeeded in producing only 
one further weak detection of the pulsar. PSRs J 1857-1027 
and J1055-6905 appear to null on more typical timescales 
of ~l-50 pulses, as shown in Figure [ll| Future observations 
of the candidate nulling pulsars from our sample will enable 
analysis similar to that of Ritchings (1976), which in turn 
will help decide whether the observed flux density variations 
arise due to nulling or simply as a result of interstellar scin- 
tillation. 



4 CONCLUSIONS 

We have conducted a survey of intermediate latitudes of 
the southern Galaxy for pulsars at ~1400 MHz. The new 
13-beam 21-cm receiver of the Parkes radio telescope was 
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used to rapidly cover to moderate depth a large region of 
sky flanking the area of the deeper ongoing Galactic plane 
survey (Lyne et al. 2000; Camilo et al. 2000). The inter- 
ference environment was formidable, however development 
of a comprehensive scheme for the rejection of pulsar can- 
didates arising from interference enabled the realisation of 
the full expected survey sensitivity of approximately 0.5 mjy 
for slow and most millisecond pulsars. The survey was highly 
successful, detecting 170 pulsars of which 69 were previously 
unknown, in a relatively short observing campaign. The new 
discoveries are not significantly more distant than the previ- 
ously known population in this region of sky, indicating that 
the success of the survey is attributable to its sampling of 
a different portion of the broad distribution of pulsar spec- 
tral indices. The detected sample, in combination with those 
of the Galactic plane and high-latitude surveys (when com- 
plete), will prove invaluable for population modelling due to 
the use of a common observing system to cover a large area 
of sky at high radio frequency. 

Among the most interesting new objects are two recy- 
cled pulsars with massive white dwarf companions (Ed- 
wards & Bailes 2001b), four with probable low-mass He 
dwarf companions, two isolated millisecond pulsars, and one 
'slow' pulsar with a very narrow pulse profile and small pe- 
riod derivative, suggestive of recycling in a scenario similar 
to those of the known double neutron star systems (Ed- 
wards & Bailes 2001a). As expected from the large Galactic 
z-height of much of the survey volume, the detected popula- 
tion of slow pulsars was relatively old and as such exhibited a 

fi-artinn nf pnlgarg slinwino; milling anrl snh-piilap mnrl- 



high : 

ulation 1 Two pulsars show very regular drifting sub-pulses 



and are analysed in detail elsewhere (Ord et al. in prepara- 
tion). 
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